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THE SURVIVAL OF LARGE DIMENSIONAL THRESHOLD 
CONTACT PROCESSES 

By Thomas Mountford-'^ and Roberto H. Schonmann^ 

Ecole Polytechnique Federale and University of California at Los Angeles 

We study the threshold 6 contact process on Z'' with infection 
parameter A. We show that the critical point Ac, defined as the thresh- 
old for survival starting from every site occupied, vanishes as d — > oo. 
This implies that the threshold 6 voter model on Z'' has a nondegen- 
erate extremal invariant measure, when d is large. 

1. Introduction and results. In this paper we study the threshold 9 con- 
tact process on the d-dimensional cubic lattice Z'^. [See, e.g., Liggett (1985) 
and Liggett (1999) for background on the area of interacting particle sys- 
tems, to which this process belongs.] The threshold 6 contact process on Z'^ 
is defined as the Markov process on {0, 1}^'' with flip rates at site x G Z'^ at 
time t>0 given by: 

• 1 flips to at rate 1. 

• flips to 1 at rate A in case there are at least 9 nearest neighbors of x in 
state 1 at time t; and at rate otherwise. 

The parameter A > is called the infection rate. The state of the process 
at each site at each time is called the spin at that site at that time. A spin 
is interpreted as a vacant or healthy site, while a spin 1 is interpreted as 
an occupied or infected site. It is well known [see, e.g.. Chapter I of Liggett 
(1985)] that such rates deflne a unique Markov process. Note that the flip 
rates above are attractive [see Chapter III of Liggett (1985)], a property 
that has many consequences, and will be used extensively in this paper. 

When 9 = 1, the threshold contact process is easier to analyze, among 
other reasons because it has an additive dual process [see Chapter III of 
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Liggett (1985) for additive duality]. (In this dual process infected sites be- 
come healthy at rate 1, and they infect simultaneously all their neighbors 
at rate A.) Also, only when 6 = 1 can every site become infected when the 
process is started from a finite number of infected sites. The behavior of 
the process with 6* = 1 is therefore not expected to be qualitatively different 
from that of the much studied (linear) contact process [see Chapter VI of 
Liggett (1985) and Part I of Liggett (1999)]. In contrast, when 6 >2, not 
only is the analysis of the model more challenging, but also the behavior 
is expected to be qualitatively distinct from that of the threshold 1 or the 
linear contact process. We defer to Section 5, item (2), a discussion of this 
point, and we refer the reader to Toom (1974), Durrett and Gray (1990), 
Bramson and Gray (1991), Chen (1992, 1994) and Pontes and Schonmann 
(2008) for background on threshold contact processes with 6 >2 and re- 
lated models. Here we will extend to every value of a result that is already 
available in the case 6 = 1. 

Threshold contact processes can also be seen as natural generalizations 
of the bootstrap percolation processes [see, e.g., Adler and Lev (2003), for 
a review], with bootstrap percolation corresponding to the case with A = oo 
(or in other words, with the down flips suppressed). 

Por coupling purposes, it is convenient to construct the process using a 
system of Poisson marks. Por this purpose, associate to each site in Z*^ two 
independent Poisson processes: one with rate 1, and one with rate A. Mark 
the arrival times of the former with symbols D (for "down" ) and those of the 
latter with symbols U (for "up"). Make these Poisson processes independent 
from site to site. Use the marks now in the obvious way, to define the process: 
A spin 1 at site x flips to when it encounters a D mark there; a spin 
at site X flips to 1 when it encounters an U mark there and at least 6 
neighbors of x have spin 1 at that time. The probability space on which 
these Poisson processes are defined will be large enough to accommodate 
the process started from arbitrary initial configurations. The same system 
of Poisson marks will be used to build various comparison processes that 
will be introduced through the paper. 

We will denote by {rj^a g x t)t>o — {'nt)t>o the process started from a ran- 
dom configuration picked according to law /x at time 0. When /x is concen- 
trated on a single configuration ^, we will use the notation (??t )t>o for the 
process. When /x is product measure with density p we will use the notation 
(^f)t>0- Similar shorthand notation will also be used in other places below. 

The point mass on the configuration with all sites in state i G {0, 1} will 
be denoted 5j. The distribution 5q is trivially invariant for the threshold 
6 contact process, when 6 > 0. By attractivity, rj^ ^ as i — > oo, where 
=^> denotes convergence in distribution, and i?, is called the upper invariant 
measure. 
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We say that the process started from the distribution fi dies out when 
r]^ ^ 5o, as t ^ oo. When this happens for every /i, we simply say that 
the process dies out. Attr activity imphes that the process dies out precisely 
when = Jq. When the process does not die out, we will say that it survives. 

Set = ¥{rj^{0) = 1). The critical point for survival is defined as 

Ac = sup{A : z/ = 60} = sup{A : ^ as t — > 00}. 

It is a standard, easy fact that Ac > 0, and it is also known that for each 
d>2 and 9 < d, Ac < 00. Durrett and Gray (1990) have an explicit proof 
in the case 9 = 2, using contour methods. But a well-known argument, that 
we present in Section 5, item (5), allows one to use the results of Toom 
(1974) and Bramson and Gray (1991) in discrete time to obtain this result 
in continuous time, for arbitrary 9 < d. Here we will prove the following. 

Theorem 1.1. For any 9, the threshold 9 contact process has 

lim Ac = 0. 

ci— +00 

In case 9 = 1 this result is not new. As pointed out in the display in page 
191 in Cox and Durrett (1991), in this case it is known that Ac < C /d, for 
an appropriate constant C . The argument given there comes from Holley 
and Liggett (1978), who used it in the case of the linear contact process. 

From the analysis of a mean-field version of the threshold 9 contact pro- 
cess, in Pontes and Schonmann (2008), it is natural to conjecture that for 
each 9 there is a constant $0 G (0, 00) such that dX^ — > as d — > 00. Our 
methods unfortunately are not sharp enough to prove this strengthening of 
Theorem 1.1, or even to obtain good estimates on how fast Ac vanishes, as 
d 00. 

Theorem 1.1 has an implication for the threshold 9 voter model on Z*^, 
defined by the following rates. At x ^Z'^ at time t>0: 

• 1 flips to at rate 1 in case there are at least 9 nearest neighbors of x in 
state at time t; and at rate otherwise. 

• flips to 1 at rate A in case there are at least 9 nearest neighbors of x in 
state 1 at time t; and at rate otherwise. 

When 9 > 1, the threshold 9 voter model obviously has 6q and 5i as two 
of its extremal invariant measures. In case 9 > 2, it is also easy to find 
infinitely many other extremal invariant measures concentrated on a single 
configuration that is a trap for the process. On the other hand, from the 
work of Cox and Durrett (1991), restated as Proposition 2.11 in Part II of 
Liggett (1999), we know that a sufficient condition for this process to have 
a nondegenerate (i.e., not concentrated on a single configuration) extremal 
invariant measure is that the threshold 9 contact process survives when 
A = 1. In combination with Theorem 1.1 above, we obtain 
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Theorem 1.2. For any 6, the threshold 9 voter model on Z"^ has a 
nondegenerate extremal invariant measure when d is large. 

In the case 6 = 1 a much more detailed result is available, from Liggett 
(1994): the threshold = 1 voter model on Z"^ has a nondegenerate extremal 
invariant measure if an only if d > 2. 

The proof of Theorem 1.1 will be split into Sections 2, 3 and 4. The proof 
will be presented in the case 6 = 2, for notational simplicity, and so that the 
reader can more easily visualize the arguments. In Section 5 we will explain 
why the proof applies to the other values of 9. Section 5 will also contain 
other results and comments related to Theorem 1.1 and its proof. 

2. Approach to the proof of Theorem 1.1. As mentioned in the Introduction, 
in this section and the following two, we suppose that 9 = 2. 

We will compare the threshold 2 contact process on with the same 
processes on subsets of Z"^, with free boundary conditions. For this purpose 
we will use the following notation. For each x G Z"^, let J\fx be the set of its 2d 
nearest neighbors, and given R C Z'^, set Mr^x = n R. Given Rc Rc Z*^, 
and C € {0, 1}^, we will denote by r/^.^ the threshold 2 contact process on 
R, started from the configuration C restricted to R. In this process, at site 
X £ R, at time t > 0: 

• 1 flips to at rate 1. 

• flips to 1 at rate A in case there are at least 2 sites in Mr^x in state 1 at 
time t; and at rate otherwise. 

Deflne G = Z^ x {0, l}*^"^. We wiU denote by Ac(G) the threshold for survival 
(started from all sites occupied) for the threshold 2 contact process on G. 
We will prove 

Proposition 2.1. 

lim Ac(G) =0. 

d— >oo 

By attractivity, Ac(Z'^) < Ac(G), so this result implies the case 9 = 2 of 
Theorem 1.1. In this section we will explain how the proof of Proposition 
2.1 will be reduced to the proof of two lemmas. 

Fix a value of A > 0. We want to prove that for d large enough, the 
process {TiQ-t)t>o survives. For this purpose, we will partition 7? x {0, 1}''"^ 
into copies of the hypercube H = {0,1}'^, indexed by the elements of I?. 
The elements of this partition will be denoted by 



Hi = {2ii,2n + 1} X {2i2,2z2 + 1} X {0, lY'^ , i = (ii.is) G Z^- 
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We will observe the process {r]G-t)t>o at times that are multiples of (P. At 
each one of these times, we will decide whether each hypercube Hi, i S Z^, 
is in a good configuration or not according to a technical definition provided 
later in this section. (For large d, the completely empty configuration will 
not be a good configuration, while the completely full configuration will be 
a good configuration.) This will be done in a way that will allow us to use 
a classic result of Toom (1974), to show that the hypercube H = Hq will be 
in a good configuration at arbitrarily large times, implying therefore that 
iVG;t)t>o survives. 

We recall now the result from Toom (1974), that will be used. Let X^^^, 
i E Z^, n G Z, be a collection of i.i.d. random variables, with P{Xi^n = 0) = e, 
P[Xi^n = 1) = 1 ~ where e € [0, 1] is a noise parameter. Define now the 
following family of random transformations, Te-„, on {0,1}^^, indexed by 
n £ {0, 1, 2, 3, . . .}. T£;o is the identity and recursively, for n = 1, 2, . . . , 

r Xi^ri, if {Ts;n-i^)i = 1, or 
[ 0, otherwise, 

for ^ G {0,1}^^ i G Z^. Denote by 1 the element of {0,1}^^ identically 1. 
Toom (1974) proved the following result [see also Bramson and Gray (1991), 
for an alternative proof] : There exists £c> such that 

(2.1) iore<£c, liininf P((r^;„l)o = 1) > 0. 

Returning to the threshold 2 contact process on G, we now explain how 
good configurations are defined. For this purpose, consider the threshold 2 
contact process on H, started from C G {Ojl}*^. We want to compare this 
process with the system of independent spins on H, that starts from the 
same configuration, and uses the same system of up and down marks, for 
its flips. (So that flips down occur at rate 1 and flips up occur at rate A, 
independently of the state of the other spins.) Call this independent flip 
process {T^H-t)t>o- Clearly, r?^.^ < vr^.^, for all t > 0. Roughly speaking, a 
good configuration ( on H, is one for which the processes (^£-.Jo<t<d2 and 

i^H-t)o<t<d'^ are likely to be close to each other. The technical notion of 
closeness that will be used is defined next. 

Definition {Thin sets). Given R,Sc Z^, and L > 0, say that 5 is L- 
thin in iZ, if 

|5nA/}j,^| <L, 

for all X £ R. 



6 



T. MOUNTFORD AND R. H. SCHONMANN 



Definition {Good configurations). Given Rc Rc Z'^, and C G {0, l}"^, 
we define the discrepancy sets by 

A^., = {xGi?:7?^.,(x)/4.,(x)} 

= {xGi?:7/^.,(x)=0,4.,(x) = l}, t>0, 

A^= U 4;*- 

t6[0,d2] 

And given also L > 0, and 6 G [0, 1] , we say that the configuration C is (L, S)- 
good on R, if 

P(A^ is L-thin on i?) > 1 - 5. 

The next lemma will be the key to providing us with good configurations. 
First we need one more technical definition. 

Definition (Almost product distributions). Given it! C Z*^, M > and 
a G [0, 1], say that a probability distribution u over {0, 1}^ is an M-almost 
a-product distribution, if it can be constructed as follows. There is a prob- 
ability space on which there is a random configuration ^ G {0, 1}^, with law 
given by an a-product distribution (i.e., each site of R is independently oc- 
cupied in ^, with probability a). On this probability space there is also a 
random set S C R that is almost surely M-thin on R. The law can then be 
obtained as the law of the random configuration defined as being identical 
to ^ on i?\ 5 and identically on S. (Note that in this definition, the set S 
does not need to be independent of ^.) 

Convention on notation and terminology. In the three defini- 
tions above, when R = H, we may omit it. For instance, "C is good" has the 
same meaning as "C is good on H," = A^, etc. 

Lemma 2.1 (Almost product distributions concentrate on good configu- 
rations). For any A > and any a £ (0, 1] there is L = L(A, a) = L(\/ (1 -|- 
A) A a) such that for any M > 0, 6 > 0, and e > 0, there is a finite do such 
that if u has an M-almost a-product distribution over {0,1}^, then 

uiiC-.C is {L,S)-good})> 1- e ifd>do. 

The quantity A/(l -|- A), that appears in the statement of Lemma 2.1, 
is the equilibrium density of the independent flip process. We will use the 
notation 



L(A) = L(A/(1 + A)) = L(A, a) for a > A/(l + A). 
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Note that in particular, for L = L{\) and arbitrary (5 > 0, the configuration 
with all sites occupied in H is (L,5)-good, if d is large. 

The proof of Lemma 2.1 will be given in Section 3. The next lemma, 
whose proof builds on Lemma 2.1, and will be postponed to Section 4, will 
allow us to use (2.1) to prove Proposition 2.1. Given L and 5, define the 
rescaling operator IZi^s '■ {0, 1}*^ {0, 1}^ as follows: 

cv (/-w - / 1' if C is (L,5)-good on Hi, 
{l<-LA^}U-<y^^ otherwise, 

i G Z^. Below >z means stochastically larger than. 

Lemma 2.2 (Rescaling). For any A > and e > there is 5 > and 
finite do such that for L = L{X), d > do, and any ( G {0, 1}*^, 

Iterating this lemma [using the Markov property of {VG-t)t>o]j refer- 
ring to the remark after Lemma 2.1, we learn that for any e > there is 
6 > such that for large enough d, 

(2.2) TZLAvh-M^) ^ T,,nl, n=l,2,3,.... 

Taking e < Sc in (2.2) and combining this estimate with (2.1), proves that 
for any A > 0, when d is large enough, the process {^Q.t)t>o is {L,6)-good 
on H at arbitrarily large times (multiples of d^), provided that L and 6 
have been properly chosen. Since the empty configuration on H is clearly 
not (L, 5)-good, when d is large enough, we conclude then that ir]G-t)t>o 
survives. This shows that Proposition 2.1 will have been proved once we 
prove Lemmas 2.1 and 2.2. 

3. Almost product distributions concentrate on good configurations (proof 
of Lemma 2.1). In order to prove Lemma 2.1, we introduce the following 
definition. A site x £ H is said to be (C, -f^)-dangerous if in the indepen- 
dent flip process (vr^.Jt>o there is some t G [0, d^] when x has fewer than K 
occupied neighbors. Define now 

= {x £ H -.x is (C, i(') -dangerous}. 

Note that the definition of D'^ refers only to the independent flip process 
(7r^.j)t>o, and not to the threshold 2 contact process (??^.j)t>o- The next 
lemma makes nevertheless the connection with that latter process. 



Lemma 3.1 (Only dangerous sites can cause trouble). For any ^ G {0, 1}^, 
any L > 0, and any K > L + 2, on the event that Dp^ is L-thin, we have 
C D'^. In particular, if is L-thin, then also A^ is L-thin. 
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Proof. We argue by contradiction. (Informally: we consider the first 
violation of the claim, and show that it cannot be a violation of the claim.) 
Consider 



(inf = oo). If the lemma were false, we would have s < d? and there would 
be some x ^ H \ D'^ with the following properties. At {x, s) there would be 
an up mark at which there would be a flip from to 1 in the process (tt^^.J, 

while in the process (t?^^.^) the spin at x would be immediately before and 
immediately after time s. This implies that immediately before s the site x 
had fewer than 2 neighbors that were occupied in the process {ffn t)- But, 
by the definition of s, for t < s, any site y that has VH tiv) '^H-tiv) ™ust 

be (^, i^)-dangerous. Therefore, if we are on the event that D'j^ is L-thin, we 
learn that immediately before s the site x had at most L neighbors where 
7r£-.j differed from and hence had fewer than L + 2 neighbors that were 

occupied in the process (vr^.J- This implies that x G D^; a contradiction. 
□ 

Note that the next lemma refers only to the independent flip process 
(7r^.Jt>o, and does not mention the threshold 2 contact process (^£'.f)t>o- 
The arguments used to prove this lemma have some affinities with ideas of 
Balogh, Bollobas and Morris (2007). 

Lemma 3.2 (Controlling dangerous sites). For any A > and any a £ 
(0, 1], there is L = L{X, a) > 0, such that for every M > 0, K > and e > 0, 
there is a finite = c?o(A, a, M, K, e) such that ifv is an M -almost a-product 
distribution over {0,1}^, then 



Proof. We will use the following notation for sites in if = {0, l}*^. For 
each Ac {1, . . . , d}, ca will denote the indicator of the set A, that is, the ele- 
ment of H that is 1 on ^ and on {1, . . . , d} \ ^. Note that the neighborhood 
of the origin in H is Nufi = {e{i}, • ■ • , ej^}}. 

By the symmetries of H , with no loss in generality, 



s = inf{tG [{),d^]:/^i<^D^^] 



dv{Q¥{D\ is L-thin) >l-e 



ifd>d'^ 




diy{C)¥{Dj^ is not L-thin) < 2'^ / di/(C)F(|D^ n A/W,o| > L) 
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Since v is an M-almost a-product distribution, 

1 dz/(C)P(e|i}, . . . , e|i} G < ^ d^«(OP(e{i}, . . . , e|i} G ^i+Af ), 

(3.2) 

where product distribution over {0, 1}^ with density a. 

For i = 1, . . . , L, consider now the set Ni = {e^ijj : j = L + 1, . . . ,d}. Note 
that these L sets are disjoint from each other, that they have cardinaUty 
d — L and that for each i, Ni C MH,eiiy ■ Let 

i;«'^+^ = inf|te[o,(i2]: X: 4,(e{..})<^ + ^4' 

i = 1, . . . , L. Since the event {e{j} € -D^-,,.^} is included in the event {7".^'^+*^ < 
d^}, and the sets Ni are disjoint from each other, 

J d/x,(OP(e{i},...,e{^| GZ)i+^,) < J d//«(e)P (^fl {^/™ < ^'}) 
(3.3) ^ 

i=l 

We want to estimate from above the probabihty that < d^. Set 

W^ = {t£ [0,(^2 + 1] :Ei=L+i4;t(e{i,i}) <K + M}, and let be the 

Lebesgue measure of this set. Define also T/ as the instant of the first up mark 
in Ni after Note that T/ — 7"^'^+^^ exponentially distributed with 

rate {d — L)X, and is independent of 7"?'^+^^, Therefore, for (i> 1/A, using 
Markov's inequality, 

.i.P(Tf-',.>.p(-(lz^) 

eXd * 

Hence, using Fubini's theorem, 

'd/.„(e)p(T;«'™'<d') 



(3.4) < eAd J d/u„(OE(|Wf I) 



< eAd((i^ + l)P(Bin(d-L,aA ) <i^ + M 
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where Bin(n,p) stands for a binomial random variable corresponding to n 
independent attempts, each with probability p of success. And we used the 
fact that in the independent flip process (vr|^.t), each site is independently 
occupied at time t with a probability that varies monotonically in t, starting 
from a at t = 0, and converging to the equilibrium value A/(l + A) as t — > 
oo (as well as the fact that binomial random variables are stochastically 
monotone increasing in the probability of success). 

In order to choose L{X, a), we recall the standard large deviation estimate 

(3.5) pfBinf^,a A — ^) < ^fa A < exp{--fd), 



where 7 = 7(0 A j^) > 0. For a reason that soon will become clear, we take 
L = L{X, a) = L(A/(1 + A) A a) sufficiently large so that 



(3.6) 

From (3.1)-(3.5) above. 



exp(-7L) < i. 



diy{C)¥{D)^ is not L-thin) 



e\d{d^ + 1) 

X p(Bm(d-L,aA 



1 + A 



<K + M 



(3.7) 



<2'^d^ 



eXd{d^ + 1) 



xP(Bm(|,«A^ 



df A 

< T r 

~ 4V 1 + A 



1 L 



< 2''d'^[e\d{d' + l)exp(-7d)]^ 
= {e\d\d^ + l)t{2eM-lL))\ 

where in the second inequality, we suppose that d > di , for some large enough 
di=di{L{\a),K + M). 

Given e > 0, it is clear that (3.6) and (3.7) imply the existence of a finite 
d'^ = (io(A, Q, M, i^, e) (larger than di and 1/A), so that the claim stated in 
the lemma holds. □ 



Proof of Lemma 2.1. Take L = L(A,a) from Lemma 3.2, and apply 
Lemma 3.2 with K = L + 2 and e replaced by eS, to obtain the existence of 
do = do(A, a, M, 5, e) such that 



j dv{C)^{D^j^ is not L-thin) < e5, 
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for d > dQ. From Lemma 3.1, we have for every ( G {0, 1}^, 

P(A'^ is not L-thin) < F{D'}^ is not L-thin). 

On the set {C : C is not {L, 5)-good}, the probabihty on the left-hand side of 
the last display is larger than 6, by the definition of {L,6)-good configura- 
tions. Therefore, integrating that inequality over u, and using the previous 
one, we obtain 

5z^({C:C is not {L,6)-good}) < J dz^(C)P(A^ is not L-thin) < £(5. 
Canceling 6, we obtain the desired claim. □ 

4. Rescaling (proof of Lemma 2.2). Set 

H = HU H(^Q i) U H(i Q) . 
The main step in the proof of Lemma 2.2 will be the following lemma. 

Lemma 4.1 (Local rescaling estimate). For any A > and e > there 

is 6 > and finite do such that for d>do the following holds. If C ^ {0, 1}^, 
satisfies one of the following two conditions: 

(a) ( is {L{X),6)-good on H, 

or 

(b) C is {L{X),6)-good on -ff(o,i) '^''T'd on H[ifi), 
then 

P(r/1^2 is {L{\),5)-good onH)>l-e. 

Lemmas 4.3 and 4.4 below will contain the basic ingredients in the proof 
of Lemma 4.1. 

Lemma 4.2 (Goodness is monotone). If C <C" ^i"^ iwo configurations 
in {0, 1}^, then C A^ . In particular, if for some L > and 6 > 0, C is 
{L,6)-good, then also Q" is {L,6)-good. 

Proof. If x G A^", then there is some up mark at x at some time 
s < d^, with the following properties. At (x,s) there is a flip from to 1 in 
the process (vr^-.^), while in the process (??^.J the spin at a; is immediately 
before and immediately after time s. This implies that immediately before 
s the site x had fewer than 2 neighbors that were occupied in the process 
(yn-t)- Since, by attractivity, ry^.^ < for all t > 0, it follows that also 
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in the process {jYn-i) immediately before s the site x was vacant and had 
fewer than 2 neighbors that were occupied. Therefore, x wiU still be vacant 
in the process (^/^.J, immediately after s, while in the process (vr^.^) it will 
be occupied. This implies that x G A'^' . □ 

Lemma 4.3 (Goodness is likely to be preserved). For any A > 0, Af > 0, 

and 5 > 0, for each ( £ {0,1}^ that is {M,6)-good and t G [d"^ /4:,d'^], we 
have: 

(i) The total variation distance between the law of rf^.^ and a certain 
M -almost \/{\ + \) -product distribution is bounded above by 25, for large d. 

(ii) For any 5' > 0, 

is {L{\),6')-good) > 1 - 45 for large d. 

Proof. Let be a random configuration in {0,1}^, corresponding to 
a product distribution with density A/(l + A). This is the equilibrium dis- 
tribution for the independent flip process {iTH,t), so that = 7r|^j has the 
same A/(l + A)-product distribution for each t > 0. Define now 

eA^:Cc \_i^t{x), if X G i7 \ or is not M-thin, 
~ \ 0, if X G A^ and A[ is M-thin. 

Note that, for arbitrary i > 0, S^f'^''^ has an M-almost A/(l + A)-product 
distribution. 

The event A that there is at least one Poisson mark, either up or down, 
at each site x £ H during the time interval [0,(i^/4] has 

P(^^) <2'^exp(-(iV4) <<5, 

for large d. And on A, vr^.^ = ,^ti for t>d'^/4. 

Let B denote the event that A^ is M-thin. From the assumption on 
we have 

P(S^) < 5. 

On AB we have 

^^^=\0 = r,Ux), if.GA^ \=Vh-A-)- 



This means that ry^-.^ = ^f^'^, on the event AB, which has ¥{{ABy) < 26, 
when d is large, completing the proof of (i). 

The claim (ii) follows from (i) and Lemma 2.1, by possibly taking d larger. 

□ 
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Lemma 4.4 (Infecting with goodness). For any A > 0, M > 0, and /3 > 0, 

if Tp is a random configuration on {0, 1}^ that has restrictions to -f^(o,i) o'^'^ 
^{1,0) that are independent and each has an M-almost jS-product distribu- 
tion, then we have: 

(i) Tj^ ^ has a restriction to H that is stochastically larger than an 8M- 

almost a-product distribution, where a = a{X,P) > 0. 

(ii) For any e > and 5 > 0, 

F{r]^ is {L{X,a),6)-good on H) >l — e for large d. 



Proof. For each i = (^1,^2) G partition Hi into the sets 
Hi^y = {2ii,2ii + 1} X {2^2,2^2 + 1} X {y}, y € {0,1} 



d~2 



From the hypothesis, we have, for i = (0, 1) and i = (1, 0), a set Si C Hi that 
is M-thin in Hi, and such that the restriction of tp to Hi is identical on 
Hi \ Si to & random configuration ^ on -ff(o,i) U -^(1,0) that has a /3-product 
distribution. 
Define 

1= {ye {0, i}'^-^ : //(o 1)^^ n 5(0,1) / or i7(i,o),, n S(i,o) / 0}, 

so that for y G {0, l}'^^^ \ I, tp \s identical to ^ on -ff{o,i),y U -H^(i,o),j/- 

For y = (2/3, . . . G {0, l}'^"^, consider the following event Ay-, ^(x) = 1 
for X S -f^{o,i),y U H(i Q-^ y, between times and 1, there is no down mark in 
H{o,o),y U H{o,i),y U -f^(i,o),y ^11^ between times and 1 there is an up mark 
at each one of the sites (1, 1, 2/3, . . . , y^), (1,0,^3,...,^^), (0, 1, 2/3, . . . , and 
(0, 0, 2/3, ... , yd), in that order, or in the order (1, 1,1/3,..., Vd), (0, 1, 2/3, ... , yd), 
(1, 0, 2/3, ■ • • , yd), (0, 0, 2/3, • • • , yd) - Then the events Ay, y e {0,lY~'^ , are i.i.d., 
with 

P(^) > P^e-^\l - e-^l^f = a'{X, (3) > 0. 

Define also a random configuration 1^ on by setting 

if 2; S H^Q Q-^ y for some 2/ for which Ay occurs, 
\ 0, otherwise. 

Set 

U^{o,o),y 

We have now the following estimate, which is central to our argument, and 



iiiitiuiin duuvc, aiiu. uiic way yi^^ 

structed from the Poisson marks: 



follows immediately from the definitions above, and the way {i]^ ) is con- 



r]t^{x)>(l){x) iix£H\S. 
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Claim (i) in the lemma will therefore follow, once we argue that the law of (f) 
is stochastically larger than a product distribution with some density a > 0, 
and that S is 8M-thin. 

For the first of these tasks, it is enough to take a = a(A,/3) = a'((A,/3)/4. 
Indeed, a product distribution over H with density a gives, for each y G 
{0,1}'^"^, probability (1 — a'/4)^ > 1 — a' to the event that all 4 sites in 
H(o,o),y ^-re vacant. And this is enough, since in (j), when H^^fl-^ y is not fully 
vacant, it is fully occupied. 

We now argue that S is an 8M-thin subset of H. For this, let and 
S"!^ be, respectively, the translates of ^(o,!) and <S'(i^o) to H. Then S' = 
5|q U S"!^ is 2M-thin. To obtain S from S' , we introduce the mappings 
ak-.H^H, A: = (A:i,A;2)G{0,l}2,givenby 

, / _ J 1 — Xj, if j G {1, 2} and kj = 1, 

- I if j G {3, . . . , 4 or kj = 0. 

These 4 mappings are automorphisms of H (as a graph), and 

S= U akiS')- 

So S is the union of 4 sets that are 2M-thin, and hence it is 8M-thin. 
This completes the proof of claim (i). 

Claim (ii) is now a consequence of (i) and Lemmas 2.1 and 4.2. □ 

Proof of Lemma 4.1. In case condition (a) holds, the claim in the 
lemma is immediate from (ii) in Lemma 4.3 [with M = L{X), 5' = 6], provided 
that we take 5 < e/4. 

In case condition (b) holds, we will show next that the claim in the lemma 
holds, provided that we take 5 < e/12. For this purpose, we will compare 
the process (r/^ )t£\od'^]j with the process (r?^ )te\o d'^] iii which interaction 
is only allowed among the hypercubes H, -ff(o,i) and -f^(i,o) during the time 
interval [(iV2- l,dV2]- 

The random configuration fi^ „ has restrictions to H, Hrn-i) and 

//(I 0) that are independent. From (i) in Lemma 4.3 [in its translated versions 
for -ff(o,i) and H(^i q^, with M = L{X)], we learn that, for large d, this random 
configuration has restrictions to -ff(o,i) and -f^(i,o) each with a law that is at 
total variation distance at most 26 from a L(A)-almost A/(l + A)-product 
distribution. Therefore the random configuration fi'x „ has a law that is 

^ 'H;d^/2-l 

at total variation distance at most 45 from a random configuration on H 
that satisfies the conditions of Lemma 4.4, with M = L{X) and (3 = A/(l + A). 
[See, e.g.. Lemma (6.2), Section 2.6, page 139 of Durrett (1996), for the fact 
that the total variation distance is added when measures are multiplied.] 



>F{7]t is (L(A,a),<5)-good on H) - 4S > 1 - e/3 - 46, 
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Using now the Markov property of {ffi^^^t&[i),(P] (apphed at time — 1), 

and (ii) in Lemma 4.4 (with the time interval [0,1] translated to — 
l,(i^/2]), we learn that for an appropriate a > 0, 

P(4.rf2/2 is {L{\,a),5ygood on H) 

for large enough d. 

Using once more the Markov property of {ff-. )ti^\Qd?] (this time applied 
at time and using (ii) in Lemma 4.3 [this time with M = L(A,a:), 

t = (i^/2 and 5' = 5], we learn now that 

P(7?1^2 is not (Z(A),(5)-good on H) 

<F(4,rf2/2 is not (L(A,a),5)-good on H) 

+ ¥{ff^ 2 is not (L(A), 5)-good on H\ff^ ^ is (L(A, a), 5)-good on H) 

< e/3 + 4(5 + 4(5 < e, 
for large enough d. The claim in the lemma now follows by using attractivity 

4.) 



and Lemma 4.2, when comparing the process {ff^ ^t&[Q,d?] with the process 

('?|.Pte [0,(^2]- n 

Proof of Lemma 2.2. We will combine Lemma 4.1, with Theorem 
O.O(i) of Liggett, Schonmann and Stacey (1997). For this purpose, for i E 1?, 
set 

Hi = Hi U i^j+(o,i) U i^j+(i^o) 

and define the random variables X^, i G Z^, by setting X^ = 1 in case 

{TIlAv^^.J^ = 1> or (7^L,^(C))^ = (^L,5(C))i+{o,i)(^L,<5(C))*+(i,o) = 0; and 

X^ = 0, otherwise. Then [X^)i^-^2 has a finite range of dependence. There- 
fore, Theorem O.O(i) of Liggett, Schonmann and Stacey (1997) assures us 
that given e > 0, there exists e > such that if 

(4.1) infP(X^ = l)> 1-e, 

then there is an i.i.d. random field {Xi)i^j2, with P{Xi = 0) = e, P{Xi = 
1) = 1 — e, that lies stochastically below (-^f )igz2- 

On the other hand, Lemma 4.1 (with i in place of the e there), assures 
us that (4.1) is indeed satisfied, provided that we take 5 > Q properly, and 
d large enough. 
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Define the fohowing random transformation, T, on {0, 1}^ : 

{f^)i = I ^i"' if = 1 or ^'i+(o,i) = ^i+{i,o) = 1, 
1 0, otherwise, 

for ^' G {0, 1}^^, i £l? . Under the conditions above on 5 and d, we have 
now, 

-R-L^rfcd^)^ TTZlAO^ Te-iTlLAO- 
This completes the proof. □ 

5. Extensions, related results, notes. We will use the following standard 
notation, for the elements of the canonical base of Z'^: ei = (1, 0, 0, ... , 0, 0), 
62 = (0,1,0,... ,0,0), . . . , erf = (0,0,0,..., 0,1). For x G Z'^, we say that its 
nearest neighbors y and z are colinear with x if {y, z} = {x — ej,x + ej} for 
some j S {1, . . . , d}. 

We will be considering other processes on Z'^ which can be constructed 
with the same system of Poisson marks used for the threshold contact pro- 
cesses. As usual, we say that model A dominates model B, if when they are 
both constructed in this way and are started from the same configuration, 
model A is at any time above model B. 

(1) First we explain how the proof of Theorem 1.1 should be modified 
for arbitrary 6. One uses the following generalization of (2.1). Let 
i G Z'^, n G Z, be a collection of i.i.d. random variables, with P{Xi^n = 0) = e, 
P{Xi n = I) = 1 — £, where e G [0, 1] is a noise parameter. Define now the 
following family of random transformations, T^-n, on {0,1}^'^, indexed by 
n G {0, 1, 2, 3, . . .}. Te;o is the identity and recursively, for n = 1, 2, . . . , 



Xi^n, if {Ts-n-l'^)i = 1, or {Ts-n-l^)i+ej = 1, j = 1, . • . , 

0, otherwise, 



for ^ G {0, 1}^ , i G Z'^. Denote by 1 the element of {0, 1}^ identically 1. 
Toom (1974) proved the following result [see also Bramson and Gray (1991), 
for an alternative proof] : There exists £c> such that 

(5.1) foTe<ec, liminfP((r,;„l)o = l)>0. 

For our purposes, we set 

G = Z''x{0,l}^-^= U Hi, 

where 

Hi = {2ii, 2ii + 1} X • • • X {2ie,2ie + 1} x {0, l}'^-^ i = {h, ■ ■ ■ ,ie) & ■ 



THRESHOLD CONTACT PROCESSES 



17 



One then proves the analogue of Proposition 2.1, exactly in the same way as 
that proposition was proved, except for obvious changes [e.g., (5.1) is used 
in its ^-dimensional version, in Lemma 3.1 the condition reads K > L -\- 6, 
H = HU {[jj^i QHej), in (i) in Lemma 4.4, the quantity 8M becomes 
02^M, etc.]. 

(2) Even when the threshold 9 contact process survives started from 
all sites occupied, it may happen that for small p> 0, rjf ^ 6o, as t — > oo. 
By attractivity, if this happens for some value of p, it will also happen for 
smaller values of p. This leads to the following definitions: 

Pc{^) = sup{p G [0, 1] : ?7f =^ 5o, as t — > oo} 

= sup{p G [0, 1] : pf* ^ as t — > co}, 

A = inf {A : PciX) = 0} = inf I A : for all p > 0, inf pf > o|. 

In case 6 = 1, one can use duality to easily show that whenever the process 
survives, it survives starting from any product distribution with positive 
density. In particular, Ac = A. 

Chen (1994) proved a result that implies that if = 2, d>3, then 

(5.2) A<oo. 

Specifically, Chen (1994) studied a spin flip model on Z*^ that he referred to 
as the system with symmetric sexual reproduction. That model is the 6 = 2 
case of the model that we call the modified threshold 6 contact process, and 
is defined by the following rates. At a; G Z'', at time t > 0: 

• 1 flips to at rate 1. 

• flips to 1 at rate A in case there are at least 6 nearest neighbors of x, 
no two of which are colinear with x, in state 1 at time t; and at rate 
otherwise. 

Theorem 4 of Chen (1994) shows that when d>3, the modified threshold 2 
contact process on Z'^, with a large enough A, survives for any initial density 
p> 0. By domination, the same therefore holds for the threshold 2 contact 
process on Z'^, that is, (5.2) holds. 

In Chen (1992) the following result, that contrasts with (5.2), was proved. 
In d = 2 the modified threshold 2 contact process with arbitrarily large A, 
started from a product measure with small enough density p> (how small 
depends on A) dies out. This means that for this model A = oo. Presumably 
the same holds for the threshold 2 contact process in dimension 2. 

Interestingly enough, the qualitative behavior of the threshold 2 contact 
process in high dimensions is distinct from the behavior of the associated 
mean- field model, analyzed in Pontes and Schonmann (2008). In that mean- 
field model, Pc(A) > for all A > 0, with \imx^oo PcW = 0. In particular, A = 
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oo. Also the behavior of the threshold 2 contact process on a homogeneous 
tree is distinct from that of both the process on Z"' with large d and the 
mean-field model. For the models on the homogeneous trees with > 2, 
Pontes and Schonmann (2008) observed that even limA^ooPc(A) > 0. 

In connection with the facts reviewed above it seems very interesting to 
decide whether, for 6 = 2, in d>3 the critical points Ac and A are distinct 
from each other. 

(3) It is worth pointing out that Theorem 1.1 extends also to the mod- 
ified threshold 6 contact process, as defined in (2) above. This is so simply 
because the modified threshold 9 contact process dominates the threshold 
26 — 1 contact process. 

(4) The following strengthening of Theorem 1.1 can be proved with very 
little additional work. 

Theorem 5.1. For any 6, and any A > 0, 

lim Pc(A) = 0. 

d— >oo 

To prove this theorem, first we prove that for arbitrary p> 0, e' > and 
6'>0, for L = L(A), 

(5.3) nL,s'irfa.a^)h T,,.^l 

for large d. This claim is a direct consequence of Lemma 2.1 [with a = 
p, M irrelevant, e = e' /2, L = L{X,p) and 6 = e' /8], combined with (ii) in 
Lemma 4.3 [with M = L{X,p), 6 = e'/8 and t = d'^], and with the fact that 
the right-hand side of (5.3) has a product distribution with density e' . 

Combining (5.3) with n — 1 iterates of Lemma 2.2 [using the Markov 
property of (??G' t)t>o]) we learn that for any e > there is 5 > such that 
for L = L{X) and large enough d, 

(5.4) nL,s{vl,^d^) t T,,nl, n= 1,2,3,.... 

Theorem 5.1 now follows from (5.4) in the same way as Theorem 1.1 followed 
from (2.2). 

(5) We end this section with a proof that must be well known, but does 
not seem to be in print. We observe how (5.1) can be used to prove that for 
every d>2, the asymmetric model defined by the rates below has Ac < oo. 
At X e U^, at time t > 0: 

• 1 fiips to at rate 1. 

• fiips to 1 at rate A in case x + ej, j = 1,. . . ,d, are all in state 1 at time 
t; and at rate otherwise. 
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We will use {fjt)t>o to denote this process. Define the random variables 
Xi^n) i £ ^'^5 n = 1,2,..., by Xi^^ = 1 if at i, during the time interval 
[{n-l)/VX,n/VX] there is no down mark and there is an up mark; „ = 0, 
otherwise. Then the Xi^n are i.i.d., and e = P(Xj^„ = 0) can be made arbi- 
trarily small, by taking A large. It is also clear that 

't/Vx^ 1,2,3,.... 

Therefore (5.1) implies survival of {f]t)t>o when A is large enough. By dom- 
ination, we learn also that for the threshold contact processes on Z'^, with 
9 < d, Ac < oo, as mentioned in the Introduction. 
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